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Abstract 


Three ongoing microlensing experiments have found more candidate events 
than expected from the known stars. These experiments measure only one param¬ 
eter of the massive compact halo objects (machos), the magnihcation time scale of 
the events. More information is required to understand the nature of the machos. 
A satellite experiment has been proposed to measure their projected transverse 
speed V = vjil — z), where v is the macho transverse speed and 2 ; its distance 
divided by the distance of the source. Measurement of v would determine whether 
the machos were in the Galactic disk, Galactic halo, or in the Large Magellanic 
Cloud (LMC). We simulate events observed toward the LMC by the Earth and 
by a satellite in an Earth like heliocentric orbit. To leading order, such an ex¬ 
periment determines h up to a two fold degeneracy. More precise measurements 
break the degeneracy. We show that with photometric precisions of 3% to 4% and 
approximately 1 observation per day, v can be measured with a maximum error 
of 20% for 70% to 90% of events similar to the ones reported by the EROS and 
MACHO collaborations. The projected transverse velocity is known with the same 
maximum error for 60% to 75% of these events. This 20% maximum error is not 
a 1 (7 error but is mostly due to degeneracy between two possible solutions, each 
one being localized to much better than 20%. These results are obtained with an 
Earth-satellite separation of 1 AU, and are improved by a larger separation. 

Subject Headings: astrometry - dark matter - gravitational lensing - Magel¬ 
lanic Clouds 
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1. Introduction 


The EROS collaboration (Aubourg et al. 1993; Aubourg et al. 1995a, 1995b) 
and the MACHO collaboration (Alcock et al. 1993, 1995b) have reported to date a 
total of hve candidate microlensing events of stellar sources in the Large Magellanic 
Cloud (LMC), more than predicted from the known luminous stars of the Milky 
Way (Bahcall et al. 1994; Gould, Bahcall & Flynn 1995) and the LMC (Gould 
1995b), but less than expected with a standard spherical halo (Paczyhski 1986; 
Griest 1991). The MACHO collaboration (Alcock et al. 1995a; Bennett et al. 
1995) and OGLE collaboration (Udalski et al. 1994) have reported several dozen 
candidate events toward the Galactic bulge, also more than expected from known 
stars (Paczyhski et al. 1994; Han & Gould 1995; Zhao, Spergel, & Rich 1995). 

To understand both the nature of the massive compact halo objects (machos) 
which generate these events and the structure of the Galaxy, one must collect as 
much information as possible about each event. The current microlensing searches 
yield only one characteristic of the lens, the time scale R of the event: R = 
where v is the magnitude of the transverse velocity v of the macho relative to the 
Earth-source line of sight, and is the Einstein radius of the macho: 

4(j 

rl = —dsourceMz{l - z). (1.1) 


Here M is the macho mass, dgource is the distance between the Earth and the 
source (in e.g. the LMC) and 2 ; is the ratio of the macho distance di^^s to the 
source distance: 


2; = 


di 


ens 


d 


source 


( 1 . 2 ) 


A variety of methods of obtaining additional information have been discussed. 
The most efficient and promising one seems to be the observation of the same 
events by both Earth-based and space-based telescopes (Refsdal 1966; Gould 1994, 
1995a). This method uses the fact that the event looks substantially different as 
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seen from the two locations to measure two new characteristics of the macho, the 
two components of the “projected velocity” v, the transverse velocity projected 
onto the observer plane, 


V 




= V- 


^ler 


V 

1 - Z 


( 1 . 3 ) 


Table 1 gives the value of the projected speed for the different possible populations 
of machos on lines of sight toward the LMC (disk, thick disk. Galactic halo, LMC 
halo, LMC disk). It shows that by measuring the projected speed of each event 
with an error of ^ 20%, one can determine the population to which each macho 
belongs with good conhdence. Thus, for example, one could determine whether the 
machos make up a signihcant fraction of the dark halo. Moreover, the additional 
measured parameters constrain the kinematics and the mass spectrum of whatever 
population is detected. 


TABLE 1 

Mean Projected Velocities 


Component (n) 2 ; (h) 


km s ^ km s ^ 


Galactic Disk 

50 

0.01 

50 

Galactic Thick Disk 

100 

0.04 

105 

Galactic Halo 

220 

0.20 

275 

LMC Halo 

100 

0.94 

1700 

LMC Disk 

30 

0.99 

3000 


By comparing the peak magnihcations and peak times of the event as seen from 
the Earth and satellite, one could determine v up to a four fold degeneracy and v up 
to a two fold degeneracy (Refsdal 1966; Gould 1994). Gould (1995a) showed that 
these degeneracies can in principle be broken by measuring the small difference in 
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the event time scales due to the relative motion of the Earth and the satellite. He 
also gave a rough estimate of the precision required to break the degeneracy. Here 
we report on Monte Carlo simulations of observations of the same microlensing 
events from the ground and a satellite. We have determined the conditions needed 
to break the degeneracy for a large fraction of the events and which would permit 
measurement of the projected speed v with an error ^ 20%. We focus attention 
on the photometric precisions, the frequency of the observations, and the Earth- 
satellite separation. The simulations consider a broad range of possible macho 
masses. We restrict our study to observations toward the EMC. Observations 
toward the Galactic bulge are physically and geometrically different, and will be 
discussed elsewhere (Gaudi & Gould 1995). 

2. Error Analysis 

The mathematical analysis underlying the Monte Carlo is substantially more 
involved than is typically the case. We therefore include in this section a complete 
mathematical description. The reader who is interested primarily in the results can 
skip directly to § 3. However, as we discuss at the end of this section, most of these 
results can be understood analytically and therefore study of the present section 
is useful for a physical as well as a mathematical understanding of the problem. 

The magnihcation A of a microlensed star is a function only of the separation x 
of the lens and star in units of the angular Einstein radius 9^ = r^jdi^ns- Explicitly 
A{x) = + 2)jx{x"^ + 4)^A (Paczyhski 1986). Assuming rectilinear motion of 

the observer, source, and lens, x evolves with time as 

x{t) = - hf + /3% (2.1) 

where ui = {3 is the impact parameter in units of the Einstein radius, and % 

is the time of maximum magnihcation. The hux F[t) from the microlensed star is 
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therefore a function of five parameters = (to, Fq, i?), 


F{t] to, /3, Fq, B) = FQA[x{t] to, /3, ^)] + i?, 


( 2 . 2 ) 


where Fq is the unlensed flux and B is the light from any unresolved additional 
sources that are not lensed (such as the lens itself or a binary companion to the 
lensed star). If a set of measurements yj^ are made at times then the parameters 
may be ht by minimizing = '^^.[yk ~ Fitk] where Uk is the error in the 

A;th measurement. By differentiating in the neighborhood of the solution one 
hnds Cjj, the covariance matrix of the Oj, 


= 


j, _ _-29F{tk-, ai...a5) dF{tk-, ai.-.a^) 

k ■' 


(2.3) 


If measurements of the event are also made from a satellite, then the flux 
can be described by hve additional parameters a) = jS'^ Fq^ B'). (In the 
simulations described in the following section, we assume that the errors are given 
by CTk = iJEFo[A{tk)]^F and cr^ = asF^lA'{tk)]^F where and cry are hxed 
parameters.) The covariance matrix for both sets of measurements combined is 
then a 10 X 10 block diagonal matrix with the upper block given by equation (2.3) 
and the lower block given by a similar equation for the satellite. There are now 10 
parameters a^. 

If the Earth and the satellite telescopes have the same hlters, then To = Fq. 
Formally, this constraint can be written 

= 0, (2.4) 

where = (0, 0, 0,1, 0, 0, 0, 0, —1, 0) and where we have introduced the Einstein 
summation convention. To see how to impose the constraint, assume hrst that it 
has a hnite but very small error Q: = 0 F Q. Then the associated inverse 
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covariance matrix is 6^ = aia^/Q'^. The restricted covariance matrix is then 
see e.g. Gould 1989). Letting Q ^ 0 yields, 


Cij — Cjj 

CmnCtmCtn 


(2.5) 


If the Liters are the same, the unlensed light is also equal, B = B' ^ and this leads 
to a similar constraint with = (0, 0, 0, 0,1, 0, 0, 0, 0, —1). Equation (2.5) should 
then be applied twice successively, once for each constraint. 

To a first approximation, the Earth and the satellite are at rest with respect to 
each other, so that io ~ u>'. As discussed by Gould (1994), the projected velocity 
is then given by 

/\~y 

v = Ax = (cuAf, A/3). (2.6) 

Here L is the magnitude of L, the projection of the Earth-satellite separation 
vector onto the plane of the sky. At = Tq — to, and A/3 is the difference in impact 
parameters. The components of Ax are respectively parallel and perpendicular to 
the projected separation vector. There is a two fold degeneracy in the magnitude 
|A/3| = A/3± depending on whether the source is on the opposite or same side 
of the macho as seen from the Earth and satellite. There is a further two fold 
degeneracy in the sign of A/3 depending on the relative orientation of the source- 
lens separation as seen from the two observers. However, these degeneracies can 
be broken by measuring the small difference Au = to' — to due to the relative 
Earth-satellite motion. This results in a constraint 


u^^uAt + Uj_A/3 — LAio = 0, 


(2.7) 


where M|| and uj_ are the components of the velocity of the satellite relative to the 
Earth parallel and perpendicular to the projected satellite-Earth separation vector 
(Gould 1995a). 
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When constructing the simulations, we of course know which side of the macho 
the source passes as seen from the Earth and the satellite and hence we know not 
only the magnitudes of fi and (3'^ but also the magnitude and sign of A/3. However, 
we initially assume that we measure only the magnitudes of /3 and (3' and that all 
four possible values of A/3 are equally likely. We then test these possible solutions 
separately. Consider for example the +A/3_ solution, which we label “—[-”• We 
hrst form the 3-vector of parameters a~~^ = (At,+A/3_, Acu). We then calculate 
the 3x3 covariance matrix of these quantities, c”"*" = Cj' ~ ~ C'-i-5,j + C+Sj+S- 

Next we formulate the constraint as in equation (2.4) with L, —1) 

which leads to an adjustment of the covariance matrix c —)• c [see eq. (2.5)] 
and an adjustment of a~~^ (see e.g. Gould 1989), 


—t _ ~—h/r—1-„—1-\ _ „—h 


= c • 




■(& 


jk 


= a- — 


aittf 


-+ 


C-mn 




( 2 . 8 ) 


where b 
constraint: 




We evaluate the of this best (—h) solution given the 


rl+ = (ar+ 




[aiaj 


-+\2 


C-mn CtmCtn 


(2.9) 


At this point in the analysis, we assume that the best estimates of At, Acu, /3 
and (3' are equal to the true values. (We consider the effect of possible differences 
between the best estimates and the true values below, but as we eventually show, 
these are relatively unimportant.) For the true values and for the right choice of 
A/3, we have ajaj = 0. Since the four possible choices of a; differ only in the second 
(A/3) component, can be written 


2 {a2S/3 +f 

X-+ = —1+ -, 

C-mn 


( 2 . 10 ) 


where 6/3 = A/3 — (+A/3_) is the difference between the true and trial values 

of A/3. For example, if the true solution is A/3 = —A/3- then h/3"*“"*“ = 2(3'^ 
8(3 = 2{f3' — /3), 6/3 = 0, and 6/3~^ = —2/3. 
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This analysis is repeated for all four solutions and is evaluated for each. The 
only important modihcation is that for the ±A/3_|_ solutions, the A/3 components of 
the covariance matrix are = ±(— C2,i + £2,1+5 — £7,2 + £7,2+5) (z ^ 2) and 

C 22 ^ = £ 2,2 + £ 2,7 + £ 7,2 + £ 7 , 7 , while for the — A/3_ solutions, the (z ^ 2) components 
are = cT” = — cf/"*"- 

IL LI iL 

In the simulations, we accept solutions if < 9. We also calculate the in¬ 
trinsic error in each of the allowed solutions. The fractional error in the vector 
solution is dehned as ((£v)^)^/^/£ which for small values can be approximated as 
((£Ax)^)^/^/Aa;, where we recall that Ax = (cuAf, A/3) = ( 07 , 02 ). To calculate 
this quantity, we apply the constraint (2.7) to e.g. c and form £ [see eq. (2.5)] 
and then compute 


V 


£77”*" + 2£]^2~*~ T £ 
. + («2“+ 


--1 


1/2 



( 2 . 11 ) 


The fractional error in the scalar solution is dehned as ((££)^)^/^/£ which can 
be approximated as {{SAx)‘^)^^‘^/Ax ~ (l/2)([£(Aa;)^]^)^/^/(Aa;)^. After some 
algebra, we hud 


V 




1/2 


,-+72 


+ (o 


.+)^ 


( 2 . 12 ) 


We can give some further insight into the formalism derived in this section 
by restricting consideration to the case at hand: observations toward the LMC. 
Since the LMC is very close to the ecliptic pole, the constraint can be very well 
approximated by a = (0,0,1) where 0 = 277 yr“^ is the frequency of the Earth’s 
orbit. Then equation (2.10) can be written as 


X 


2 


ilteS/3 


2 


1 + 




(2.13) 


where is the dimensionless error in the time-scale difference, ap = C 22 i 

and where we have for simplicity assumed that the correlation coefficient is small. 
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For typical parameters (/3,/3^ ~ 0-5), one finds ~ 0.5 for the ibA/3_ 

solutions and cr^ _|_/(T^ ~ 8 for the ibA/3_|_ solutions. As discussed by Gould (1995a), 
the higher errors in the latter case arise from the strong correlation in the errors 
of fi and /3'. This correlation is induced by the constraints Fq = Fq and B = B'. 

Consider hrst the ±A/3_ cases. For G ^ cTlo jjS ~ 116 days, 


X-± = 


ftteSiS 


-±\ 2 




(2.14) 


One might in principle also consider the opposite regime. First, however, such 
long events are not observed. Second, even if they were observed, they would be 
susceptible to ground-based parallax (Gould 1992; Alcock et ah 1995c) so that 
space-based parallax would not be necessary. Thus, equation (2.14) holds for all 
relevant cases. The situation is quite different for the ±A/3_|_ cases. Here 

xl±=(— —(G^Tdays), (2.15) 

V / 

so that the degeneracy breaking is independent of time scale for all but the shortest 
events. 

Using equations (2.14) and (2.15), one can gain a good idea of the sensitivity 
of the entire experiment. First, note that for hxed /3 and (3'^ and for hxed number 
of observations per Finstein radius crossing time, the errors in Ax = (cuAl, A/3) 
and in Acu/cu are also hxed. The distributions of /3 and (3' depend only on the ratio 
dsatjf'e of the Farth-satellite separation to the projected Finstein ring. Hence, if 
one ignores the problem of degeneracy breaking, the sensitivity of the experiment 
depends only on this ratio and not on the time scales of the events. The principal 
problem in degeneracy breaking is to separate the A/3_|_ from the A/3_ solutions. 
If this can be done, it is possible to determine v. If the true solution is A/3_ then 
the A/3_|_ solutions must be eliminated and visa versa. Clearly, the former is more 
difhcult because the A/3_|_ errors are larger. Thus, the overall requirements of the 
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experiment are set by equation (2.15). Notice that this is also independent of tg 
(provided tg ^ 7 days). The smaller value of min= 2min{/3,/3^} ~ (9(1). It 
can of course be less than unity, but if it is too much less, then the A/3± solutions 
become so similar that it is not important to break the degeneracy. 

Hence, one arrives at a basic requirement ^ 0.1 to routinely break the 
principal degeneracy at the 3(7 level. This implies ^ 0.01. If this requirement 
is met, then it follows immediately that the intrinsic errors around each solution are 
unimportant compared to the degeneracy breaking errors provided dgatj^e ^ 0.05. 
Thus one expects the experiment to attain a maximum sensitivity over a fairly 
broad range from dgat ^ 7g (beyond which the satellite falls outside the Einstein 
ring) and dgat ^ 7g/20 (beyond which the intrinsic errors dominate and rise as 

Finally, one can see that if the experiment is designed primarily to break the 
(two fold) speed degeneracy, then the (four fold) velocity degeneracy will usually 
also be broken. The most difficult case is ±A/3_. For dehniteness, take +A/3_. 
Then the ratio of the y^’s for the vector to scalar degeneracy breaking is 


X-- 


O'/,) 


_\ 2 


h/3++) 


tg 




22 days 


h/3++ ) 


( 2 . 16 ) 


Again, while it is not impossible that SjS be much smaller than the typical 
values of if it is too small, then the ±A/3_ solutions are so similar that it is 

unimportant to break the degeneracy. 

These qualitative conclusions are born out quantitatively in the remainder of 
the paper. 
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3. Halo Model and Monte Carlo 


The precise structure of the Milky Way halo is still unknown. In the interests 
of simplicity, we perform our simulations with the simplest existing model: the 
isothermal sphere. We discuss the effects of using other models below. 


We hx the macho mass M for each simulation and take the macho mass density 
p to be, 

do 


p{r) = 


1 + (r/ao)2’ 


(3.1) 


where r is the distance to the Galactic center (assumed to be 8 kpc at the Sun) and 
ao = 5.6 kpc is the core radius. The results do not depend on the value of po- The 
macho velocity distribution is taken to be an isotropic Gaussian in the Galactic 
frame. The distribution of the transverse speed v (in the plane perpendicular to 
the line of sight) is 

dN f \ 

— oc/(,.)= « exp (-^ . (3.2) 


where \/2a = 210 km s ^ is normalized to the local Galactic rotation speed. 

Although their effect is small, we take into account the motion of both the 
Sun and the LMC in the Galactic frame. Both are ~ 230kms“^ (Jones et ah 
1994). These motions induce a motion of the line of sight, and a small correlation 
between the probability distributions of the position and the speed of the machos. 
We neglect this correlation. The contribution to the event rate of a macho scales 
as p(r) Cg n oc p(r) ^ z{l — z) v. Hence we choose randomly the macho position in 
the line of sight with the probability distribution ^z(l — z) p(r), and the macho 
speed V with the probability distribution v f[v). The angle between the direction 
of this velocity and a hxed direction in the plane perpendicular to the line of sight 
is chosen randomly from a uniform distribution between 0 and 27r. 

The impact parameter /3 is either hxed, or chosen randomly with a hat distri¬ 
bution between 0 and fimax- We hx the Earth-satellite separation, the precisions 
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of the ground and the satellite photometries, and the frequency of observations. 
To simulate the experiment, we assume that the hrst satellite observation starts 
when the source enters the Einstein ring as seen from the Earth, and we assume 
that both sites stop observing when the source-lens separation is greater than Scg 
for both the Earth and the satellite. The events where the minimum source-lens 
separation for the satellite (3' > f.bcg are assumed to be unresolvable. 

We then calculate the four possible solutions for the projected velocity v, and 
their intrinsic errors. Eor each solution, a value quantihes its compatibility with 
the degeneracy breaking constraint [eq. (2.7)]. Those solutions with < 9 are 
retained, and those with y^ > 9 are excluded. If the fractional difference between 
all the allowed solutions and the true value is less than 20%, and if the intrinsic 
fractional errors (at 1 cr) of the allowed solutions are less than 20%, we consider 
that the degeneracy is broken and that we know the projected speed with an error 
% 20%. (As we discuss analytically in § 2 and numerically in § 4, the total errors are 
dominated by the degeneracy breaking, so the total error is still typically % 20%.) 
Eor each set of parameters, the simulation is repeated 5000 times, yielding Monte 
Carlo estimates with statistical errors less than 1%. 

The effects of using a different model can be understood qualitatively as follows. 
Eor models (such as a flattened halo model) in which the machos are a factor ( 
closer to the Earth, the effect is to make the Einstein ring smaller by This 

is very much the same as reducing the mass by a factor ( and the effect can be 
judged by examining Eigure 3 below. Eor models with slower macho velocities, 
the effect is very similar to having larger Earth-satellite relative velocities. This 
tends to increase the accuracy as can also be seen in Eigure 3 below. The principal 
model for which one would expect substantially different velocities is a truncated 
halo model, or more generally, a halo whose density falls r“", [n > 3). Eor such 
models, the velocity dispersion is lower than for an isothermal sphere by y^ 2/3 
(see problem 4.9 of Binney & Tremaine 1987). One also expects the machos to be 
closer to the Earth in such models, so there is little net effect on the degeneracy 
breaking fraction. 
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4. Results I: The Speed Degeneracy 


The simulated degeneracy breaking fraction (DBF) depends on six parameters: 
the photometric precisions cr^ and cry of the Earth and satellite measurements, the 
number of observations done per Einstein radius crossing time the maximum 
value l3max of /3, the hxed mass M of the machos, and the Earth-satellite separation 
dsat- (Because we assume the satellite is in a heliocentric Earth-like orbit, we 
measure this separation in time, one month being equivalent to an angle of 27r/12.) 
An observing rate Ugig = 20 corresponds to one observation per 0.5 to 4 days 
for the events reported to date which are typically A ~10-75 days. We therefore 
consider this to be a realistic observing rate. Fiducial parameters for the whole 
simulation are crj^ = 3%, cry = 4%, Ugig = 20, dgat = 2 months, M = O.IMq, and 
l3max = 0.7; in this case we hnd DBF = 80%. Here, we present curves representing 
the DBF when 1 or 2 out of the 6 parameters are changed. 

Figure 1 presents the DBF against Ug^g for photometric precisions cry; and cry 
of respectively 4% and 5%, 3% and 4%, and 2% and 3%. To resolve 75% of the 
events, 25 observations per Einstein ring crossing time are necessary in case of 4% 
and 5% photometry, 15 for 3% and 4% photometry, and only 8 observations are 
sufficient for 2% and 3% photometry. Hence it is possible to break the two fold 
ambiguity for a large majority of the events with realistic observing conditions. 

Figure 2 presents the DBF against /3, which is not randomly chosen but con¬ 
stant for each point of the curve. The events where /3 > 0.7 are much harder to 
resolve. They should therefore be separated in our statistical analysis. For 0 < /3 < 
0.7, the mean DBF equals 80%; however for 0.7 < /? < 1, it equals 15%. 

Figure 3 presents the DBF against the macho mass Af, which is constant for 
each point of the hgure, for hve different Earth-satellite separations. The observed 
microlensing events reported up to now correspond very probably to masses in the 
range 0.01-1 Mq. We focus on these masses in the hgure. We show curves for dgat 
= 0.5, 1, 2, 4, and 6 months because one possibility is to launch a satellite which 
would drift away from the Earth at a roughly constant rate of several months per 
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Figure 1. Degeneracy Breaking Fraction as a function of observation rate nghs (number of observations 
per Einstein radius crossing time) for three levels of fractional photometric accuracy from the Earth {(Te) 
and the satellite (cs). Shown are ((Te,<^s) = (2%,3%) {solid curve), (3%,4%) (short dashed curve), and 
(4%,5%) (long dashed curve). Other parameters are held Rxed at M = 0.1 Mq, dgat = 2 months (1 AU), 
and fdfjiax — 0.7. 


year. The set of curves then represents the evolution of the satellite efficiency over 
time. 

The two fold ambiguity in v can be most easily broken when dgat is close to 
the projected Einstein radius fe'. when dgat ^ the event looks very similar from 
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Figure 2. Degeneracy Breaking Fraction as a function of dimensionless impact parameter /3. 
parameters are held fixed at nghs = 20, (Je = 3%, as = 4%, M = 0.1 Mq, and dgat = 2 months. 


the Earth and the satellite, so the satellite observations are of little help; when 
dsat ^ f'e^ the probability that the satellite will not observe the event is high. 
Since be oc there is an optimal macho mass for each dgat, the mass for which 

the DBF is a maximum on the curves, and this optimal mass increases with dgat- 
Since it is the relative velocity of the Earth and the satellite that allows one to 
break the degeneracy, and since this velocity increases with dgat^ the maximum 
DBF also increases with dgat- 


Other 
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Figure 3. Degeneracy Breaking Fraction as a function of logio of the macho mass in solar units for 
hve values of dgat, the distance between the Earth and the satellite. Shown are dgat = 0.5 months (0.26 
AU) {dot-dashed curve), 1 month (0.52 AU) (long dashed curve), 2 months (1 AU) (short dashed curve), 4 
months (1.73 AU) (solid curve), and 6 months (2 AU) (bold solid curve). Other parameters are held hxed 
at nobs = 20, as = 3%, as = 4%, M = 0.1 Mq, and /Smax = 0.7. 


The satellite efficiency depends weakly on dgat at ~ 0.01 Mq, and increases 
quickly with dgat in the range 0.1-1 Mq. It is important to notice that one always 
benehts by increasing dgat- Globally, the DBF remains higher than 40% and most 
of the DBF values are in the range 60-90%. 
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Frequency 


For very large masses (beyond the window shown in Fig. 3), the events are 
also well resolved: for M = IOMq, the DBF is between 50% and 80% when dgat 
is larger than 1.5 months. The results are less satisfactory for very small masses: 
when M = O.OOIMq, the best DBF values are in the range 20-30%, obtained when 



0 .05 .1 .15 .2 .25 

Intrinsic Error 


Figure 4. Number of trials as a function of the largest intrinsic error in v [see eq. (2.12)] of any allowed 
solution. Trials where there were two allowed solutions which differed by more than 20% are excluded from 
this plot. Parameters are held Rxed at Ughs = 20, (Te = 3%, as = 4%, M = 0.1 Mq, [3max = 0.7, and 
dgat = 2 months. 


18 




In § 2, we argued analytically that the uncertainty in determining v is domi¬ 
nated by the problem of resolving the degeneracy between A/3± solutions and that 
the intrinsic error around each solution plays only a minor role. Figure 4 provides 
a numerical justihcation for this claim. It shows the intrinsic error for all cases 
where either a) there was only one allowed solution or b) there were two or more 
allowed solutions, but these differed by < 20%. Note that the intrinsic errors are 
typically ~ 2% and are rarely greater than 10%. 

5. Results II: The Velocity Degeneracy 

As we emphasized in the introduction, the most important requirement for 
a parallax satellite is that it break the two fold degeneracy in projected speed 
h, because this allows one to determine to which component of the Galaxy the 
detected machos belong. In the previous section we determined the observational 
requirements to meet this goal. 

Nevertheless, the additional information contained in the projected velocity v 
can also be important. For example, if the machos lie in a non-rotating isotropic 
halo, then their velocities should be characterized by a projected asymmetric drift 
~ 100 km s“^. On the other hand a rotating halo might have little or no asymmetric 
drift. 

In this section we test how well the velocities are measured given the observa¬ 
tional parameters = 20, cr^ = 3%, cry = 4%) required to break the speed 

degeneracy as determined in the previous section. For degenerate solutions, we 
dehne the velocity error as the magnitude of the difference between the true and 
the allowed solution, divided by the true speed. We evaluate the intrinsic errors of 
each solution using equation (2.11). As in the previous section, we demand that 
each type of error is less than 20%. The results are presented in Figure 5 and 
can be directly compared with Figure 3. We hud that the DBF for velocities is 
only slightly reduced from the DBF for speeds. For example, for Af = 0.1 Mq and 
dsat = 2 months, the respective fractions are 67% and 80%. In other words, it 
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Figure 5. Vector Degeneracy Breaking Fraction as a function of logio of the macho mass in solar units 
for hve values of dgat- Symbols are the same as in Fig. 3. Vector degeneracy breaking shown here requires 
that the projected velocity v be determined to 20% whereas the DBF shown in Fig. 3 requires only that 
V be determined to this accuracy. 

should be possible to extract substantial information about the directions of the 
machos as well as their speeds. 
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6. Conclusions 


We find that in feasible observing conditions, with a ground photometric preci¬ 
sion of 3% and a satellite photometric precision of 4%, and when the Earth-satellite 
separation is larger than one month (i.e., 0.5 AU) it is possible to measure the pro¬ 
jected speed V of events with C 20% accuracy for at least 70% of LMC microlensing 
events similar to the ones reported by the EROS and MACHO collaborations. 

Under the same conditions we hnd that the projected velocity v can be resolved 
at the same level at least 50% of the time. 

In these analyses we have deliberately excluded events with impact parameters 
fi > 0.7. Such events are not at present being discovered by EROS and MACHO 
because they are essentially excluded by the event-selection criteria. If such events 
are discovered in the future, they should be excluded from satellite follow-up ob¬ 
servations because, as we have shown, the probability that the reduced speed can 
be measured is very small (15%). 

We hnd that it is generally advantageous to have larger rather than smaller 
Earth-satellite separations: there is greater precision for masses M ~ 0.1-1 Mq 
and there is nearly no penalty for masses M ~ 0.01 Mq. Since it is expected that 
large separations are also favored for observations toward the bulge (Gould 1995a; 
Gaudi & Gould 1995), large separations appear to be favored by all considerations. 
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